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Abstract 

We compare a number of methods used to locate resonances. These include the speed plot, the 
time-delay method of Eisenbud and Wigner, the time-delay matrix of Smith, and a modification 
proposed by Ohmura. Numerical results show a consistency not previously reported. One time- 
delay method gives the most consistent results, supporting the conclusions of a recent theoretical 
study. 

PACS numbers: PACS numbers: 11.55.Bq, 11.80.Et, 11.80.Gw 



1 



Given a set of elastic scattering amplitudes or, for a multi-channel process, the full S- 
matrix, the most direct way to locate resonances involves a search for poles in the complex 
energy plane. To do this, one must have an analytic representation that can be continued, 
correctly accounting for the branch cuts. If only numerical values of the amplitudes are 
available, less sophisticated methods may still applied, at least for those resonances with a 
clear resonance signature. 

Here we will compare two such methods — the speed-plot, advocated by Hohler and 
the time-delay, which was first studied for elastic and multi-channel scattering by Eisenbud 
and Wigner j^, [J. These methods have recently been examined in terms of their theoretical 
basis 0,0]. Here we will briefly review the proposed methods and make comparisons based 
on their ability to correlate speed or time-delay peaks with resonance energies. In doing 
so, we will gain a better understanding of some puzzling results associated with previous 
applications of these methods. One particular time-delay method is advocated in the recent 
study of Ref . [5] . A comparison of our numerical results appears to support their conclusions. 

In 1948, Eisenbud [2J proposed a method for locating resonances, in both single- and 
multi-channel scattering, based on the time-delay associated with a wavepacket. Wigner 
quoted 3| a time-delay result in a later paper focused on the restrictions of causality. The 
result of Eisenbud was 

Ai E = h^[arg(S-l% (1) 
which for a single channel, gives the result 

^ = n%, (2 , 

S being the phase shift. The result quoted by Wigner was larger by a factor of 2. A reason 



for this factor was later noted by Wigner in Ref. [6, 

Smith [8j derived a time-delay matrix, based on the flux passing through an interaction 
region of radius R. His result for the average lifetime of a metastable state due to a collision 
beginning in the i th channel was 

Q = ~^% s ^ (3) 

S being the S-matrix including all open channels. Smith then claimed that his result could 
be connected to the Eisenbud result, using the following representation 



Atij S = Re 



-ih(Si, 



-l dS{j 

He 



(4) 
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implicitly attributed to Eisenbud. From the above two results, Smith noted that 

Q = EW ! . (5) 



followed trivially. Notice, however, that Eqs. (1) and (4) are not equivalent. 

Objections to Smith's time-delay formalism were expressed in the extensive study of 
Ohmura 9] as noted in Ref. 6J. Ohmura's result agreed with that of Eisenbud, though 
this was not apparent as the Eisenbud thesis remained unpublished and the Eisenbud result 
was misquoted by Smith. Ohmura further claimed that Smith's results could be made more 
consistent by subtracting off a term due to the outgoing unscattered wave, giving for the 
average time delay [9j of a wavepacket beginning in the i th channel 



< Ati >av = Nu 1 Re 



- ih J2(S in * - 6 in )^£ 
ah, 



with 



Ohmura's result can be found if one starts with the expression 

* d , , AT dT* 
^{argT} = T- - T-, 

uses Eq. 1, and multiplies by appropriate factors of H and 2i, to obtain 



Re 



~ ikSij * dE 



(6) 



(?) 



(8) 



(9) 



for i ^ j. This can be compared to Smith's result in Eq. 5. The sum differs in the % — j 
term, which is constructed to be consistent with Eq. 1, in the single-channel scattering limit, 
rather than Eq. 4. However, this sum is not a proper average, since the weights do not add 
to unity. To have the sum represent the average of Eisenbud time-delays, the weights must 
be renormalized, and the factor Nu in Eq. 7 is added for this purpose. 

In the following numerical comparison, we will be using the partial-wave decomposed 
S-matrix. However, the results of Eisenbud were derived assuming a single dominant partial 
wave, whereas the Ohmura result was based on the use of the full S-matrix . It is thus 
surprising how well these methods work, applied to a set of partial wave amplitudes derived 
from a fit to data. 

In Figs. 1 and 2, we compare these time-delay methods, and the more common speed-plot 
method, applied to amplitudes obtained in a fit to elastic pion-nucleon scattering. Figure 1 



compares the time delay results of Smith and Ohmura, based on the scattering S-matrix. In 
the fit to pion nucleon elastic scattering and r/N production data of Ref. loj], a multi-channel 
Chew-Mandelstam K-matrix formalism was used. This produced a multi-channel S-matrix, 
with each channel having the required poles and cuts (for r)N, 7rA and pN). As the 7rA and 
pN channels were constrained only by the ttN channel inelasticity, the displayed results are 
a test of the method and would be improved using a more detailed multi-channel analysis. 

Figure 1(a) shows the Sn channel, containing the N(1535) and N(1650) resonances, plus 
a cusp at the r]N threshold. Peaking in the time delay is evident at both resonances, 
denoted by arrows at the real parts of the pole positions U|. Separation of the N(1535) 
and the threshold cusp is clear using both forms. In Figure 1(b), results for the Pn (Roper) 
channel are given. Finally, in Figure 1(c), the D 13 resonance is shown. In these plots, the 
normalization factor Nu is proportional to the imaginary part of the irN elastic T-matrix. 
For the Roper resonance in particular, dropping this factor causes a significant shift in the 
peak position, as the imaginary part of T n N is rapidly increasing above and below the energy 
corresponding to the real part of the pole position. 

In Figure 2, we display the speed plot, which is given by the absolute value of dT/dW, 
and compare this to the single-channel result of Eisenbud, which is given by the energy 
derivative of the phase of the T-matrix. Again, for orientation, we locate the real parts of 
the associated pole positions. Peaks in the speed plots again correspond to the pole positions 
and the Eisenbud time-delay peaks. Separation of the rjN cusp and N(1535) resonance is still 
evident, though not as clear as in Figure 1. (The region of negative time-delay between the 
resonances is not plotted.) In the study of Hohler, the 7]N threshold and N(1535) resonance 
signal were combined into a single peak at the r]N threshold. The separation is visible in 
our case only if a very fine grid of energies is used. 

All methods agree when applied to the elastic P33 partial wave and the A(1232) reso- 
nance. All methods find a peak corresponding to the pole position. We have also applied 
these methods to the higher partial waves and generally find peaks for all PDG 4-star reso- 
nances using either the S-matrix methods of Smith and Ohmura or the single-channel results 
obtained usin g th e Eisenbud time-delay or speed plot. These results are quite different from 
those of Ref. 12], finding no prominent peaks for most of the isospin 3/2 resonances. In 



that work Eq. 4 was used for both elastic and inelastic resonances. 

Problems with resonance location and the use of Eq. 4 are clearly illustrated if we examine 
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the 4-star A (1950) resonance, which produces a nearly canonical resonance loop in the 
Argand plane. This resonance has T^n/Tt < 1/2 and therefore the loop passes below the 
center of the Argand circle. The four methods applied in Figs. 1 and 2 are compared to 
Eq. 4 in Fig. 3. All produce peaks at the same point, the real part of the energy associated 
with the pole, except Eq. 4, which produces a sharp dip closer to resonance mass found in 
a Breit-Wigner fit to the amplitude. In Figs. 4 and 5 of Ref. [ijj], resonance positions were 
associated with the positive time-delay shoulder, which is also evident in Fig. 3, prior to the 
time-advance spike. 

One 4-star resonance, the A(1620) occurring in the S31 partial wave, breaks the pattern 
described above. The amplitude, time-delay, and speed plots for this partial wave are 
displayed in Fig. 4. In an Argand diagram, this amplitude first moves in a clockwise direction 
before beginning the (anti-clockwise) resonance loop. Both matrix methods and the speed 
plot produce peaks near the resonance energy. However, the Ohmura result is shifted from 
the pole position. Both of the single-channel time-delay results fail to peak at the resonance. 

In summary, we have examined a number of candidate methods for the location of reso- 
nances using either the T-matrix for a single channel, or the S-matrix accounting for all open 
channels. These methods generally produce a peak at an energy associated with the real 
part of the pole position and show remarkable agreement in most cases. One method which 
has been used extensively, based on the energy derivative of the channel phase shift, does 
not correlate with the pole position, producing a peak at the pole mass of the A (1232), but 
time-advance dips near the Breit-Wigner masses of some other more inelastic resonances. 
The A(1620) resonance appears to provide the most demanding test (from a well established 
state). For this state, the matrix methods and speed plot peak near the resonance energy, 
the single-channel time-delay methods do not. Apart from the speed plot, only the S-matrix 
method of Smith appears to consistently correlate with resonance pole positions. This agrees 
with a recent claim from Kelkar and Nowakowski Q| who have reviewed the theoretical basis 
for the methods used in the present numerical study. We expect these result will motivate 
further investigations parallel to those reported in Ref. Q]. 
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FIG. 1: Time-delay plots based on the S-matrix approaches of Smith 8] (dashed) and Ohmura 
(solid) in arbitrary units. Arrows denote the real part of resonance pole positions for (a) the Sn, 
(b) Pn, and (c) D13 partial waves (see text). 
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FIG. 2: Speed plots (solid) and Eisenbud single-channel time-delays (dashed) in arbitrary units. 
Arrows denote real parts of resonance pole positions for (a) the Sn, (b) Pn, and (c) D13 partial 
waves (see text). 
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FIG. 3: Time-delay and speed plots for the F37 ttN resonance A(1950) in arbitrary units. Ar- 
rows denote the real part of resonance pole position. Plotted are (a) the real(solid) and imag- 
inary(dashed) parts of the dimensionless F37 ttN amplitude, (b) the matrix average time delay 
of Smith (dashed) and Ohmura (solid), the result using Eq. 4 (dot-dashed), (c) the Eisenbud 
single-channel time delay (dashed) and the speed (solid). 
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FIG. 4: Time-delay and speed plots for the S31 irN resonance A(1620) in arbitrary units. Ar- 
rows denote the real part of resonance pole position. Plotted are (a) the real(solid) and imag- 
inary(dashed) parts of the dimensionless S31 7riV amplitude, (b) the matrix average time delay 
of Smith (dashed) and Ohmura (solid), the result using Eq. 4 (dot-dashed), (c) the Eisenbud 
single-channel time delay (dashed) and the speed (solid). 
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